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Abstract 

We extended the Broad Histogram Method in order to obtain spectral degeneracies for 
systems with multiparametric Hamiltonians. As examples we obtained the critical lines for 
the square lattice Ising model with nearest and next-nearest neighbor interactions and the 
antiferromagnetic Ising model in an external field. For each system, the entire critical line is 
obtained using data from a single computer run. We also discuss the accuracy and efficiency 
of our method. 



PACS numbers: 02.70.Lq, 05.50.+q, 75.10.Hk 

The search for new and more efficient methods for computer simulations is always an intensive 
task in science. Multispin coding techniques (see |jj and references therein), cluster algorithms 
0, [| f|], reweighting procedures |5], |6], [7| and methods which calculate directly the spectral de- 
generacy m |], [10], 1 1] are a few remarkable examples (see [12, 13, [14|] for reviews). An efficient 



implementation of these methods can lead to an enormous increase in speed and accuracy of com- 
puter simulations. Multiparametric simulation is another way to further improve the efficiency. 
Multiparametric simulations allow, from a single computer run, to explore the whole space of 
parameters defining a given system, instead of repeating the whole process every time some pa- 
rameter is changed. 

For single-parameter simulations, several methods such as the histogram method || , simulated 
tempering ]7| , multicanonical ensemble || ||, 1/fc-sampling f|lQ| , broad histogram [JTTJ] among 



others have been shown to be successful. On the other hand, only a few works in the literature treat 
multiparametric simulations using such methods. The histogram method || is easily generalized 
to any number of parameters of the Hamiltonian. However the main problem in this approach is 
that the number of simulations needed to cover a given region in the space of parameters increases 



1 



with the system size, in a <i-dimensional space, as L d l 2 for each parameter ||. Shteto et al. used 
the multicanonical method (Entropic Sampling formulation) in order to study relaxation paths in 
the simulation of an Ising model in an external field flBf . In that work they pointed out that the 



entropic sampling can be generalized to any number of parameters. Although easily generalizable, 
a simple implementation of this method has been shown to display a very slow convergence H|. 



On the other hand, methods such as the traditional multicanonical formulation || seem to be not 
easily extendable to perform simulations with more than one parameter. Recently a re-formulation 
of the broad histogram method was introduced including an external field as a second parameter 
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In this paper we perform a spectral degeneracy calculation by applying the Broad Histogram 
method (BHM) to multiparametric Hamiltonians. We show that the method can be generalized 
to any number of parameters and the full range of them can be exploited. BHM was introduced 
three years ago [JO]] and it is based on an exact relation between the spectral degeneracy and some 



special macroscopic quantities defined within the method itself. A remarkable feature of BHM is 
its generality, since it is not restricted to any specific dynamical rule: the macroscopic quantities 
needed can be obtained by different procedures. BHM has been applied to a variety of magnetic 
systems with accurate results in a very efficient way [TTL [IB], [T7|, [TS], [Tj], |2T], |23| . 



Any statistical system can be described by two different classes of parameters: the ones which 
characterize it and are defined in the Hamiltonian, and the others which control the interaction 
with the environment. Let us consider the Hamiltonian of a given system as composed by P 
independent interaction terms 

W(X) = 7i5i(X) + 72 S 2 (X) + . . . + lP S P (X) (1) 

where the terms Si,...,Sp depends only on the microstate X of the system. 71,..., 7p are the 
parameters associated to these terms and all together define the complete Hamiltonian. For 
instance, Si could represent the coupling with an external field, S2 a first-neighbor interaction, and 
so on. The temperature, on the other hand, belongs to the class of parameters which determines 
the interaction with the environment. 

Before presenting the method, it is fundamental to introduce the concept of macrostate. To 
a given state X of the system, each <Sj(X) gives a value s». These values define a vector s = 



(si, S2, ■ ■ ■ , Sp). Hence, the set of all states X corresponding to the same s is called a macrostate. 
The space defined by all macrostates is called space of macro states. The interaction with the 
environment corresponds to the choice of the statistical ensemble, described by a weight function 
f(H, cyi, a p.) = /(71, 7p, cki, ••■) a P f i s ) whose functional form is known a priori. The control 
parameters for such interaction are a±, ap f . Once the system (characterized by s, 71, 7p) and 
its interaction with the environment (characterized by ai, Oip f ) are defined, the thermodynamic 
average of any quantity Q, can be written as 



where g(s) is the spectral degeneracy, i.e. the number of possible microstates of the system 
corresponding to the same macrostate s. (Q{s)) is the uniform average of Q in the macrostate 
s. The parameters introduced by / are irrelevant for the definition of s and consequently to the 
calculation of g(s) (this idea has already been used in another context |18|, |24|). Also, g(s) is 
independent of the particular values (71, ...,7p) of the Hamiltonian parameters. Therefore, if one 
is able to provide the spectral degeneracy g(s) and (Q(s)), all information about the observable 
Q can be obtained for the full range of the parameters (71, 7p, «i, cup.). 

In our generalization of BHM the spectral degeneracy is calculated through the steps: 

Step 1: Choice of a reversible protocol of allowed movements in the space of microstates. 
Reversible protocol means that each allowed movement X id — > X new corresponds to another also 
allowed movement X new — > X id. This definition is independent of the probabilities prescribed to 
these movements by the particular dynamics to be used in a practical computer implementation. 
These movements are virtual ones, since they are not actually performed. 

Step 2: To compute the number iV(X, As) of allowed movements, for a configuration X, that 
changes the current macrostate s to another s + As. Then (N(s, As)) is the average of iV(X, As), 
i.e., the mean number of possible movements from macrostate s to macrostate s + As; 

Step 3: Due to the reversibility condition (step 1), the total number of possible movements 
from macrostate s + As to macrostate s is equal to the total number of possible movements from 
macrostate s to macrostate s + As. Thus, we can write down the equation 



(Q) 



E s 0(s)(Q(s))/(7i, -, TP, ai, -, ap f , s) 
Es5'( s )/(7i, -,1p,®i, -,az P ,s) 



(2) 




(3) 
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This relation has been shown to be exact for any statistical model and spectrum degeneracy . 
Our multiparametric formulation remain valid based on the same arguments. Therefore, equation 
(0) can be rewritten as 

ln 9(s + A S )-ln g ( S ) = ln ( ^ s A ^ s)> . (4) 

This equation can be iteratively solved for all macrostates s, provided (N(s, As)) is known (again, 
this relation is independent of the way (N(s, As)) is obtained). This set of equations is overdeter- 
mined. However, the spectral degeneracy can be obtained without solving all equations simultane- 
ously, since all values of the projections of As along a given direction of the space of macrostates 
are equivalent. 

The Broad Histogram relation (|3]) is independent of the procedure by which (N(s, As)) is 
obtained. For the single parameter case, virtually any procedure can be adopted in this task as, 
for instance, a non-biased random walk PI , a microcanonical simulation [^(| or the random walk 
performed by the Entropic Sampling in the space of microstates |TJ|. For the multiparametric 
case, the choice of the dynamics is crucial to an efficient implementation, because the larger 
the number of parameters, the larger and sparser is the space of macrostates. In this work we 
choose a dynamics based on the Entropic Sampling method (ESM). Lima et al. |16[ have shown 



that this dynamics associated to BHM gives better results than the traditional Multicanonical 
calculation. The dynamics is implemented as follows: first, we implement the classical algorithm 
of ESM, performing the visitation in the space of macrostates. Additionally, after a given number 
of random walk steps (in our case one Monte Carlo Step (MCS)), for each visited state X, we store 
the values of iV(X, As) cumulatively into s-histograms. This dynamics is particularly interesting 
due to its simplicity. The spectral degeneracy is calculated by the BHM relation (Q), instead, 
avoiding the shortcomings of the traditional determination of the spectral degeneracy by ESM. 

The procedure described above is quite general and can be applied to any statistical model. As 
examples of applications of the method we choose two traditional multiparametric problems: i) 
The Ising model with nearest and next-nearest neighbor interactions and ii) the antiferromagnetic 
Ising model in an external field. Within the canonical ensemble, two parameters come from the 
Hamiltonian and one (the temperature T) from the ensemble definition. 

We start describing our results for the Ising model with nearest and next-nearest neighbor 



interactions in a L x L square lattice. The Hamiltonian of this system is given by 

nn nnn 

= -J X E X -J 2 E 2 (5) 

where cr, = ±1. The summations are carried over the nearest and next-nearest neighbors, re- 
spectively. Each macrostate is defined by the vector s = (Ei,E 2 ), which is independent of the 
set of parameters J\, J 2 and T. The adopted protocol of movements is the single spin-flip. In 
this case, the allowed As are: (-8,-8), (-8,-4), (-8,0), (-8,4), (-8,8), (8,-8), (8,-4), 
(8, 0), (8, 4), (8, 8). For the implementation of the random walk in the space of macrostates within 
ESM, we need to define the number of entropy updates and the number of Monte Carlo steps 
(MCS) between each update. We choose the number of MCS between successive updates as in- 
creasing linearly with the number of updates already performed, such that when the process is 
finished, each macrostate is visited on average v times. The number of updates was chosen as 
100 and v as 4000. Fig. (1) presents the entropy S(Ei,E 2 ) = hig(Ei, E 2 ) for a 10 x 10 system, 
obtained by solving (4) from the measured averages (N(s, As)). We have not used any additional 
procedure to optimize our simulations, as for example, to use the entropy obtained from another 
simulation as initial input for the Entropic Sampling-based random walk. We also choose to store 
informations for all possible values of E\ and E 2 . Once g(Ei,E 2 ) is obtained, we can calculate 
any thermodynamical quantity, without resorting again to computer simulations. Of particular 
interest is the line defined by the critical temperature T c as a function of the ratio K = J 2 j J\. 
We have calculated the internal energy, magnetization, specific heat and magnetic susceptibility 
as continuous functions of both K and T. We have used the temperature where the specific heat 
presents a maximum as an estimative to the critical point for each K. In figure (2), we present our 
results for different lattice sizes. The lines are normalized such that at K = 0, T c = 1. In the same 
figure, we present the infinite lattice results obtained by Fan & Wu using low-temperature series 



expansions [25|. Finite size effects are noticed only near to K — —0.5 (insets). An important 
advantage of our method is that it is possible to blow-up a given region in the space of parameters 
without need of additional runs. All points in the critical lines were calculated from the same 



g(Ei,E 2 ) obtained from only one simulation for each size. Traditional approaches 23] need a 
new simulational run for each pair (T,K). To our knowledge, the results presented here are the 



most detailed available. Shteto et al. |15| suggest 32 x 32 as the limit for their multiparametric 
simulations of the Ising model in an external magnetic field with Entropic Sampling. We have 
reached this limit with a much more elaborated problem. 

As a second example, we choose to study the Ising model in an external magnetic field. The 
Hamiltonian of this system is given by: 

H = ~J^2 CTiCTj + h^2oi 

nn i 

= -JE + hm (6) 

where cij = ±1. The macrostates are defined by s = (E,m). Fig. (3) presents the entropy 
S(E, m) for a 10 x 10 system. For this model, it is interesting to obtain the critical temperature 
as a function of the applied field, for the antiferromagnetic case (J < 0). In fig. (4) we present 
our results. The theoretical line obtained by Wang and Kim is also shown [^] . The agreement is 
remarkable for all values of the field. Again the pseudo-critical temperature was determined by 
the peak position in the specific heat. Due to the excellent agreement reached already with the 
sizes we have studied, we did not simulate larger lattices. 

Concerning the CPU time necessary to perform such simulations, we present in table I the 
mean visitation on each macrostate, and the total CPU time. From these results the CPU time r 
to obtain, in average, v visitations on each macrostate scales like L a where a pa 6 as expected, since 
a factor L dxP comes from the volume dependence of the space of macrostates and L d from the 
increase in the time to perform one Monte Carlo step. Spectral degeneracy calculation methods 
(as the BHM and Multicanonical) are not based on thermodynamical concepts and hence not 
subjected to critical slowing down, large energy barriers and other problems. Therefore, it is 
misleading to compare the CPU time of our method with other Monte Carlo canonical methods. 
Additionally, for traditional methods the larger the system the larger is the number of simulations 
at different temperatures one needs in order to determine the peak position for the specific heat. 
On the other hand, within the BHM once one already knows g(Ei, E 2 )(g(E,m)), it is possible 
to obtain results for any values of J±, J2, T (J, m, T) continuously. Other quantities, as the 
free energy, which are hard to obtain by traditional methods, can also be obtained from only one 
simulation. In this sense, our method introduces a huge computational speedup. 

In conclusion, we have shown that the Broad Histogram Method can be extended to obtain 
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spectral degeneracies for multiparametric Hamiltonians. This generalization allows us to calculate 
the thermodynamical averages for any set of the parameters which define the model continu- 
ously, from only one computer run. Conversely, traditional Monte Carlo simulations need to 
perform a new computer run each time any parameter of the system is changed. Compared with 
traditional Monte Carlo simulations, this feature of BHM represents an enormous computational 
gain. Moreover, the generalization is such that the method could be applied to any statistical 
model. 

This work was partially supported by CNPq, CAPES and FAPERJ (Brazilian Agencies). ARL 
acknowledge S. L. A. de Queiroz for very useful discussions. 
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FIGURE CAPTION 



Fig. 1 - Entropy for a 10 x 10 Ising model with nearest and next-nearest neighbor interactions. 
The axes represent the nearest (Ex) and next-nearest neighbor (E 2 ) terms of the corresponding 
Hamiltonian. The entropy does not cover the whole space because some pairs (Ei,E 2 ) do not 
correspond to accessible states. 

Fig. 2 - Comparison between simulations of the Ising model with nearest and next-nearest 
neighbor interactions (for L x L lattices) and low-temperature expansion [25| (solid line). The 
symbols correspond to L equal to: 10 (empty circles), 16 (pluses), 20 (crosses), 32 (diamonds). We 
also show some Metropolis results for L=32 (filled circles). Finite size effects are more remarkable 
near to K = —0.5 (upper-left inset). In the other inset (down-right) we show the results obtained 
using the ferromagnetic susceptibility to determine the critical point, instead of the specific heat. 

Fig. 3 - Entropy for a 10 x 10 Ising Model in a external magnetic field. The axes represent 
the spin coupling term (E) and the external field term (m). 

Fig. 4 - Comparison between simulation and theoretical results for the antiferromagnetic Ising 
Model in a external magnetic field. The agreement between simulation and theoretical results is 
excellent. 
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TABLE 



Model 


L 


Mean Visitation 


CPU Time (s) 


Ji, J2 


10 


4945.2 


1191.6 


Ji, J2 


16 


4537.2 


22118.2 


Ji, J2 


20 


4425.0 


86137.8 


Ji, J2 


32 


4277.4 


1848197.1 


J,h 


10 


4761.4 


1082.4 


J,h 


20 


4356.8 


71002.0 



Table 1: CPU time and mean visitation for each lattice size. The simulations were carried out on 
a Dec Alpha 500. 
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Fig. 4 - "Broad Histogram...", A.R. Lima, P.M.C. de Oliveira and T.J. P. Penna 




